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A common ploy to reduce the switching current and energy dissipation in spin-transfer-torque driven magne- 
tization switching of shape-anisotropic single-domain nanomagnets is to employ magnets with low saturation 
magnetization Ms and high shape- anisotropy. The high shape-anisotropy compensates for low Ms to keep 
the static switching error rate constant. However, this ploy increases the switching delay, its variance in the 
presence of thermal noise, and the dynamic switching error rate. Using the stochastic Landau-Lifshitz-Gilbert 
equation with a random torque emulating thermal noise, we show that pumping some excess spin-polarized 
current into the nanomagnet during switching will keep the mean switching delay and its variance constant 
as we reduce Ms, while still reducing the energy dissipation significantly. 
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I. INTRODUCTION 

Spin-transfer-torque (STT) is an electric current- 
induced magnetization switching mechanism that can ro- 
tate the magnetization axis of a nanomagnet by exerting 
a torque on it due to the passage of a spin-polarized 
currentii^— . The STT-mechanism is routinely used to 
switch the magnetization of a shape-anisotropic nano- 
magnet from one stable orientation along the easy axis to 
the other~, and has been demonstrated in numerous ex- 
periments involving both spin- valves^ and magnetic tun- 
nel junctions (MTJs)^. MTJs, consisting of an insulat- 
ing layer sandwiched between two ferromagnetic layers 
(one hard and the other soft), are becoming the staple of 
nonvolatile magnetic random access memory (MRAM)^iS 
(see Fig. [Ij. Switching the soft layer of an MTJ with 
the STT-mechanism (STT-RAM) allows for high integra- 
tion densities, but usually requires a high current density 
(> 10^ A/cm^) resulting in significant energy dissipa- 
tioi*iS. 

One way to decrease energy dissipation in STT-driven 
switching is to fashion nanomagnets out of materials with 
low saturation magnetization Ms (e.g., dilute magnetic 
semiconductors). The spin-polarized switching current 
Is, that delivers the spin-transfer-torque and switches 
the magnetization, varies as Af^ (see Refs. [l|, Illl - [l3j ). 
so that the power dissipation I^R {R = resistance of 
the nanomagnet) should vary as M^ if R does not 
change. However, reducing Ms decreases the in-plane 
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; Magnet's plane: ^ = t 90° 

FIG. 1. Simplified schematic diagram of an STT-RAM mem- 
ory element. The nanomagnets (NMs) are on the y-z plane 
and are shaped like elliptical cylinders. NM-1 is magnetically 
hardened along the z-axis so that its magnetization direction 
is fixed. The magnetization direction of NM-2 can be rotated 
with an in-plane spin polarized current that delivers a spin- 
transfer-torque. The magnetization orientation of the free 
layer NM-2 with respect to the z-axis (0° and 180°) encodes 
logic bits and 1. 



shape anisotropy barrier E[, that separates the two sta- 
ble magnetization states along the easy axis. This hap- 
pens because E^ is proportional to the product of M^ and 
the demagnetization factor of the nanomagnet, which de- 
pends on the degree of shape anisotropy. The decrease 
in Eh increases the probability of random switching be- 
tween the two stable states, which is ^exp [— E'^/ZcTp^"— 
at a temperature T. Therefore, if one reduces Ms, then 
one must also increase the in-plane shape anisotropy (or 
aspect ratio of the magnet) commensurately in order to 
keep the barrier Ei, and the static error probability un- 
changed. Increasing the shape anisotropy, or aspect ra- 
tio, has another beneficial effect; it decreases the resis- 
tance R in the path of the switching current Ig if the 
latter flows along the in-plane hard axis of the nanomag- 
net. This further reduces the power dissipation I^R. 

It therefore appears that reducing Ms, while increasing 



shape anisotropy to keep Ej, constant, is always benefi- 
cial. There is however one caveat. Reducing Ms makes a 
nanomagnet more vulnerable to thermal fluctuationsi^ 
and can increase both the thermally averaged (mean) 
switching delay and the standard deviation in the switch- 
ing delay due to thermal fluctuations. This has a dele- 
terious effect on clock speed and clock synchronization 
in memory or logic technologies utilizing spin-transfer 
torque mechanism. Consequently, memory and logic de- 
vices utilizing materials with low saturation magnetiza- 
tion often work at low temperatures, even if the Curie 
temperature of the nanomagnet exceeds room tempera- 
ture, just so that thermal agitations are suppressedi^^. 
In this paper, we show that a better approach to contend 
with thermal fluctuations, when M^ is reduced, is to still 
work at room temperature, but use slightly more switch- 
ing current than that required by the Is oc Af^ scaling 
law. This will decrease the mean switching delay and 
its variance, while still maintaining a significant energy 
saving due to the reduced Mg. 



II. MODEL 

We study the magnetization dynamics of a nanomag- 
net subjected to a spin-transfer-torque at room tempera- 
ture by employing the stochastic Landau-Lifshitz-Gilbert 
(LLC) equatio n^"'^^ . It describes the time evolution of 
the magnetization vector in the presence of spin-transfer- 
torque, the torque due to shape anisotropy, and an ad- 
ditional random torque due to thermal fluctuations. We 
choose the dimensions of the nanomagnet such that it 
has always a single ferromagnetic domain^^. Thermal ef- 
fects in magnetization dynamics have been studied both 
theoretically^"— and experimentallyS^Ti^. 

We consider a nanomagnet (see Fig. [1]) in the shape of 
an elliptical cylinder whose elliptical cross section lies in 
the y-z plane with its major axis and minor axis aligned 
along the z-direction and the y-direction, respectively. 
The dimension of the major axis is a, that of the minor 
axis is b, and the thickness is I. The magnet's volume 
is ri = {TT/4)abl. Let 9{t) be the angle subtended by 
the magnetization axis with the -|-z-axis at any instant 
of time t and (f>(t) be the angle between the -l-cc-axis and 
the projection of the magnetization axis on the x-y plane. 
Thus, d{t) is the polar angle and (j){t) is the azimuthal 
angle. Note that when </> = ± 90°, the magnetization 
vector lies in the plane of the nanomagnet. 

The potential energy of an isolated unperturbed shape- 
anisotropic single-domain nanomagnet is the uniaxial 
shape anisotropy energy given by 



demagnetization factor expressed a: 



>32 



EsHAit) = ^M^nN.it), 



(1) 



Nd{t) = Nd^,,cos^0{t) + Nd-yysm'^e{t) sin'^4>{t) 

+ Nd^^^sin^0{t)cos'^(j){t) (2) 

with Nd-zz, Nd-yy, and Nd-xx being the components of 
Nd along the z-axis, y-axis, and a;-axis, respectively. The 
expressions for these quantities can be found in Ref. 
and they are constrained by the following relation: 



Nd-ZZ + Nd-yy + N, 



= 1. 



(3) 



We have assumed that the use of a properly balanced 
synthetic antiferromagnetic fixed layer can eliminate the 
net effect of dipole coupling on the free layer—. 

At any instant of time, the total energy of the unper- 
turbed isolated nanomagnet can be expressed as 



E{t) ^ EsHA{t) - B{t)sin^6{t) + C 



(4) 



where 



B{t) = B{c^{t)) = ^ M^n [Nd-xxCos^<l>{t) 



2 



+ Nd-yySin^(f>{t) - Nd-zz] , (5) 

(6) 



The in-plane shape anisotropy energy barrier height (us- 
ing <j) — ±90°) can be expressed as 



Es 



HA,in— plane 



= (^ M^nN, 



do 



(7) 



where iV^o = [Nd-yy — Nd-zz]- Note that the in-plane 
shape anisotropy energy barrier height is independent of 
time t even though EsHj^{t) is not. 

The magnetization M(t) of the single-domain nano- 
magnet has a constant magnitude at any given temper- 
ature but a variable direction, so that we can represent 
it by the vector of unit norm nm(i) = M(i)/|M| = er 
where er is the unit vector in the radial direction in spher- 
ical coordinate system represented by (r,0,(/)). The other 
two unit vectors in the spherical coordinate system are 
denoted by e^ and e^ for 9 and 4> rotations, respectively. 
The coordinates {9,4>) completely describe the motion of 
M(t)3. 

The torque acting on the magnetization within unit 
volume due to shape anisotropy is 

TE(t) 



-n^{t)xWE[0{t),4,{t)] 
-{2B(t)sine{t)cose{t)}e^ 



{Boe{t)sin9{t)}ee, 
(8) 



where 



Mo 



Boeit) = Boem)) = ^ M''sn{Nd-xx-Nd-yy)sin{2c^{t)). 

(9) 
Passage of a constant spin-polarized current Ig through 
the nanomagnet generates a spin-transfer-torque that is 
given byS 



where Ms is the saturation magnetization and Nd is the 



Tstt(0 = ssin9(t)e0, 



(10) 



where s — {h/2e)r]Is is the spin angular momentum 
deposition per unit time and 77 is the degree of spin- 
polarization in the current Is . 

The effect of thermal fluctuations is to produce a ran- 
dom magnetic field h(t) expressed as 



h{t) = h^{t) Gx + hy{t) Gy + h:,{t) §2 



(11) 



where hx(t), hy{t), and /iz(i) are the three components 
in X-, y-, and z-direction, respectively. We will assume 
the same properties of the random field h(i) as described 
in the Ref. [35;]. Accordingly, the random thermal field 
can be expressed as^ 



h,(t) 



2akT 



|7|(l + a2)MyAt 



G 



(0,1) 



(t) {i = x,y,z) 



.(12) 
where a is the dimensionless phenomenological Gilbert 
damping constant, 7 — 2/iB/io/fi- is the gyromag- 
netic ratio for electrons and is equal to 2.21 x 10^ 
(rad.m).(A.s)~^, /is is the Bohr magneton, My = 
/ioMsJl, 1/Ai is the attempt frequency of the ran- 
dom thermal field affecting the magnetization dynamics; 
therefore At should be chosen as the simulation time-step 
used to solve the coupled LLG equations numerically, and 
G/Qi\{t) is a Gaussian distribution with zero mean and 
unit standard deviation. Note that the variance in the 
random thermal fields is inversely proportional to the 
saturation magnetization Ms] therefore a lower satura- 
tion magnetization augments the detrimental effects of 
thermal fluctuations. 

The thermal torque can be written as 

TTH(t) = Mv n^{t) X h{t) = Peit) e^ - P^{t) eg (13) 

where 

Pe{t) == Mv[ha:{t) cose{t) coscj){t) + hy{t) cos9{t)3in(j){t) 

-h,{t)sm9{t)], (14) 

P^{t) = Mv[hy{t) cos(j>{t) - h^{t) sin(j>{t)]. (15) 

The magnetization dynamics of the single-domain 
nanomagnet under the action of various torques is de- 
scribed by the stochastic Landau-Lifshitz-Gilbert (LLG) 
equation as 



dn^{t) f dn^{t) 



dt 



dt 



(l + a2)<^'(i) = -^[a(Soe(i)-s) 
+ 2B{t)cos0{t) - [sin0{t)]-^ {Pg{t) - aP^(t))]. (18) 

The application of an in-plane spin-polarized current 
Is to produce spin-transfer torque results in an energy 
dissipation I^Rt, where R is the in-plane resistance of 
the elliptical cylinder given hy R — p{2/tt) (b/a) I with p 
being is the resistivity of the material used and r is the 
switching delay. 

Furthermore, because of Gilbert damping in the nano- 
magnet, an additional energy E^ is dissipated when the 
magnetization axis in the nanomagnet switches from one 
orientation to the other along the easy axis. This energy 
is given by the expression Ed = L Pd{t)dt, where r is the 

switching delay and Pd{t) is the dissipated power given 
by22i38 



Pdit) 



a\j\ 



(1 + a^)Mv 



|TE(t)+TsTT(i)|' 



(19) 



Thermal torque does not cause any net energy dissipation 
since mean of the thermal field is zero. 

The thermal distributions of 9 and (j) in an unperturbed 
magnet are found by solving the Equations ([T7)) and ^TE\\ 
while setting Ig = 0. This will yield the distribution of 
the magnetization vector's initial orientation (Omit, 4'init) 
when stress is turned on. We consider magnetization in- 
tially situates ai 6 = 180°. The ^-distribution is Boltz- 
mann peaked at 6* = 180° with mean ~175.5°, while the 
(/)-distribution is Gaussian peaked at = ±90°—. 

The quantity r for any switching trajectory is deter- 
mined by solving the coupled equations (IT71) and p^ 
starting with an initial orientation {9init, 4'init) and ter- 
minating the trajectory when 9{t) reaches a pre-defined 
9 final, regardless of what the corresponding (j) final is. 
The time taken for a trajectory to complete (i.e., for 
9{t) to reach 9 final) is the value of r for that trajectory. 
The average value (t) and the standard deviation (At) 
are found by simulating numerous (10,000) trajectories 
in the presence of the random thermal torque, and then 
extracting these quantities from the distribution. 

The total energy Etotai dissipated in completing any 
trajectory is given by Etotai = Ed + I^Rt. The average 
power dissipated in completing any trajectory is simply 
Etotai /t. We can find the thermal average of Ed and its 
variance by calculating Ed for numerous trajectories and 
then computing these quantities from the distribution. 



hi 

Mv 



[TE(t) + TsTT(t) + TTH(i)] • (16) ,,,. SIMULATION RESULTS 



In spherical coordinate system, with constant magni- 
tude of magnetization, we get the following coupled equa- 
tions for 9- and 0-dynamics. 

(1 + a2) 9'{t) = ^[{Boe{t) - s)sin9{t) 

My 

- 2aB{t)sin9(t)cos9{t) + {aPe{t) + P4>{t))] (17) 



We consider a nanomagnet made of CoFeB alloy which 
has low saturation magnetizatioiiiS, and a low Gilbert 
damping factor of a = 0.01. The saturation magnetiza- 
tion can be varied by varying the alloy composition^^. 
The resistivity is assumed to be the same as that of 
cobalt, i.e. p = 5.81 x 10"^ f^-m [Ref. 15]. We choose this 
material over dilute magnetic semiconductors which have 
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FIG. 2. Variation of the major axis (a) and the minor axis (6) 
with the saturation magnetization (Ms) needed to maintain 
a constant in-plane shape anisotropy barrier of 0.8 eV or ~32 
kT at room temperature. The thickness of the nanomagnet 
is held constant at I = 2 nm. 
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&, and /) guarantee that the nanomagnet remains in the 
single-domain limit^iSS,. The thickness I is held constant 
at 2 nm. 

Fig. [5] shows how the major axis a and the minor 
axis b should vary with Ms to keep the in-plane shape 
anisotropy energy barrier constant at 0.8 eV. This en- 
sures that the static error probability associated with 
spontaneous switching between the two stable states 
along the easy axis remains constant as we vary M^ and 
NdQ- In Fig. 131 we plot the three components of demag- 
netization factor for different values of Mg that will keep 
the in-plane shape anisotropy barrier constant at 0.8 eV. 
Obviously, as Mg is decreased, we need to increase the 
value of Ndo, i.e. {Nd-yy — Nd-zz), to keep the same in- 
plane shape anisotropy energy barrier height. With de- 
creasing Ms, the quantity Nd-yy increases significantly 
while Nd^zz remains more or less constant. Since the 
three components of the demagnetization factor are con- 
strained by the relation Nd~xx + Nd~yy + Nd~zz = 1, the 
value of Nd-xx niust decrease proportionately, which is 
seen in Fig. [31 



A. Dependence of switching delay and energy dissipation 
on saturation magnetization 



FIG. 3. Demagnetization factors needed for different values 
of saturation magnetization Ms in order to keep the in-plane 
shape anisotropy barrier constant at 0.8 eV or ~32 kT at 
room temperature. The demagnetization factors are com- 
puted from the major and minor axes values in Fig. [2] and 
with constant thickness of 2 nm. 



much smaller Ms because the latter 's Ms is so smalli^ii^ 
that it will be impossible to make the in-plane shape 
anisotropy barrier Ei, (which is proportional to M'^Q.Ndo) 
large enough (32 kT or 0.8 eV) without making the vol- 
ume 17 of the nanomagnet very large. With that large 
volume, the nanomagnet will no longer be single-domain. 
Choosing CoFeB allows us to work at room temperature 
with a barrier of 0.8 eV or ~32 kT, while still ensuring 
single-domain behavior because the volume can be kept 
small. 

In order to maintain a constant value of E^ = 0.8 eV 
as we vary Ms , we increase the shape anisotropy of the 
nanomagnet (or the aspect ratio a/b of the ellipse) to in- 
crease Ndo and compensate for any decrease in Ms- As 
we vary the aspect ratio a/5, we keep the cross-sectional 
area of the ellipse [(7r/4) ab] and the thickness I constant, 
which keeps both the area and the volume of the nano- 
magnet n constant. The rationale behind keeping the 
cross-sectional area of the nanomagnet constant is to keep 
the density of devices per unit area on the chip constant. 

The in-plane shape anisotropy energy barrier depends 
on three quantities: Mg, Ndo and n [see Equation ([7])]. 
Since we keep J7 constant, we compensate for any de- 
crease in Ms by commensurately increasing Ndo alone. 
At all times, we ensure that the dimensions chosen (a. 



We assume that when a spin-polarized current is ap- 
plied to initiate switching, the magnetization vector 
starts out from near the south pole {6 ~ 180°) with a 
certain {9init,4'init) picked from the initial angle distri- 
butions at 300 K2i. The magnetization dynamics en- 
sures that 9 continues to rotate towards 0°, while tem- 
porary backtracking of 9 may occur due to random ther- 
mal kicks. Thermal fluctuations can introduce a spread 
in the time it takes to reach 6* ~ 0° but cannot pre- 
vent magnetization to reach ~ 0°. When 9 becomes 
< 4.5°, switching is deemed to have completed. A mod- 
erately large number (10,000) of simulations, with their 
corresponding {9init,4'init) picked from the initial angle 
distributions, are performed for each value of saturation 
magnetization to generate the simulation results in this 
subsection. The magnitude of the switching current Is is 
2 mA at Ms = 8 x 10^ A/m and it is reduced proportion- 
ately with the square of Ms for other values of Ms ■ The 
spin polarization of the current is always 80%. 

Figs, m and |5l show the mean switching delay and the 
mean energy dissipation for different values of saturation 
magnetization Ms when the in-plane shape anisotropy 
energy barrier is held constant at 0.8 eV by adjusting the 
nanomagnet 's shape (and thus demagnetization factors) 
as we vary Mg- In generating these plots, the magnitude 
of the in-plane spin-polarized current is chosen as 2 mA 
when Ms = 8 x 10^ A/m and the switching current Is was 
decreased in accordance with the Is oc M^ scaling law. 
Each scattered point (denoted 'single run') in the figure is 
the switching delay for one representative switching tra- 
jectory at that value of Ms- There is considerable scatter 
in the 'single run' data which is significantly reduced by 
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FIG. 4. Switching delay r at room temperature (300 K) in 
a nanomagnet of fixed in-plane shape anisotropy energy bar- 
rier of 0.8 eV as a function of saturation magnetization Ma. 
Switching current Is is 2 mA at Ms = 8 x 10^ A/m and it 
is reduced proportionately with the square of Ms for other 
values of Ms- The spin polarization of the current is always 
80%. 
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FIG. 5. Total energy dissipation in a nanomagnet of fixed in- 
plane shape anisotropy energy barrier of 0.8 eV as a function 
of saturation magnetization Ms at room temperature (300 K). 
The total dissipation includes the dissipation in the switching 
circuit IsRt and the internal energy dissipation Ed- The 
internal energy dissipation Ed is only a very small fraction of 
the total energy dissipation Etotai- 



thermal averaging (averaging over many trajectories) . 

The large scatter in the 'single run' data points is not 
caused by the random thermal torque since we get the 
similar trend without incorporating thermal fluctuations. 
The scatter is more prominent at smaller values of M^, 
which corresponds to lower Ig {Is oc M^). At lower Is, 
the magnetization dynamics is more complex since there 
are more ripples (see Fig. |8] and Fig. |9] later). As a result, 
there is more variability in the switching dynamics (and 
hence switching delay) with changing Is when the latter 
is small. This variability contributes to the scatter. 

In Fig. m we find that the mean switching delay (t) 
decreases with increasing saturation magnetization Ms- 
This can be explained as follows. The spin-transfer 
torque is proportional to Is/Ms- Since Is oc M^, 
the spin-transfer-torque becomes proportional to Ms- 
Therefore, reducing Ms weakens the spin-transfer-torque. 
Since the in-plane shape anisotropy energy barrier is in- 
variant, it takes longer for the weakened spin-transfer 
torque to overcome the in-plane shape anisotropy barrier 
and cause switching. This makes switching delay increase 
with decreasing Ms- 



FIG. 6. Internal energy dissipation Ed for different values of 
saturation magnetization at 300 K. This energy dissipation 
is several orders of magnitude smaller than the energy dissi- 
pation in the external circuitry thus it is only a very small 
fraction of the total energy dissipation Etotai- 



In Fig. [51 we plot the thermal means of the energy 
dissipation Etotai at 300 K as a function of the satura- 
tion magnetization Ms, while keeping the in-plane shape 
anisotropy barrier constant. Etotai is overwhelmingly 
dominated by the component I^Rt, and the internal en- 
ergy dissipation Ed has a minor contribution [see Fig.|6]. 
The switching current Is varies as the square of Ms, so 
that /g varies as M^. Furthermore, if we reduce Ms, 
we have to increase the shape anisotropy (or the aspect 
ratio a/6) to keep the in-plane shape anisotropy energy 
barrier constant. If the switching current flows along the 
minor axis of the elliptical nanomagnet (always prefer- 
able since it results in minimum resistance in the path 
of the current), then increasing the ratio a/b decreases 
the nanomagnet 's electrical resistance proportionately. 
Thus, both Is and R will decrease with decreasing Ms 
(the latter because the in-plane shape anisotropy barrier 
is kept constant). Consequently, the power dissipation 
I^R increases with Ms more rapidly than M^. Unless 
the switching delay r has a stronger dependence on Mg 
than T ex M~^, we will expect the energy dissipation to 
decrease with decreasing Mg and that is precisely what 
we observe in Fig. [51 

The last two figures highlight two important facts: (1) 
the energy dissipated to switch can be reduced by de- 
creasing Ms while maintaining a fixed in-plane shape 
anisotropy energy barrier to keep the static error prob- 
ability fixed, and (2) the switching delay increases if we 
reduce Ms while keeping the in-plane shape anisotropy 
energy barrier fixed. Thus, there are two penalties in- 
volved with reducing energy dissipation by lowering Mg 
and scaling Is quadratically with Ms'- (i) slower switch- 
ing, and (ii) higher dynamic error probability due to an 
increased variance in thermal field [see Equation ([T^ ]. 

For the nanomagnet that we have considered (with the 
parameters described earlier), we find that lowering the 
saturation magnetization Ms by a factor of ~2 decreases 
the energy dissipation by ~28 times while increasing the 
switching delay by approximately twice- The factor of ~2 
decrease in Ms causes a ^16- fold decrease in I^ (since 
Is oc M^). Additionally, there is ~4-fold decrease in 



the resistance of the nanomagnet owing to the fact that 
the shape anisotropy is increased to keep the in-plane 
shape anisotropy energy barrier constant. Thus, the ~64- 
fold decrease in power dissipation and the 2-fold increase 
in switching delay together cause a net decrease of ~28 
times in the total energy dissipation. Therefore, if we 
decrease the saturation magnetization by a factor of ~2, 
then we will: (1) gain 28-fold in energy dissipation; (2) 
lose 2-fold in switching speed; and (3) lose somewhat in 
error rates due to thermal agitation since the variance in 
switching delay is increased. 



B. Constant switching delay scaling 

In order to understand how we can maintain a constant 
switching delay while scaling Ms , let us consider the rela- 
tionship between switching current and switching delay. 
In Figs. [7] and O we plot the magnetization dynamics 
without considering any thermal fluctuations during the 
switching when M^ = 8 x lO'^ A/m and M^ = 4.09 x 10^ 
A/m, respectively. The switching current has been de- 
creased from 2 mA for M ^ = 8 x 10^ A/m to 523 ^A for 
Ms = 4.09 X 10^ A/m, in accordance with the square-law 
scaling Is oc M^. In Figs. [7]and[8l we have assumed the 
same initial orientation of magnetization Binn — 175.5° 
and (pinit = 90°, which are thermally mean values for 300 
K to avoid the stagnation point exactly along the easy 
axis. The square-law scaling however results in an in- 
creased switching delay since the latter has obviously in- 
creased by a factor of 2 (from 1.05 ns to 2.1 ns). This has 
happened because of more ripples generating from more 
precessional motion of the magnetization vector seen in 
Fig. [51 In order to maintain the same switching delay 
of 1.05 ns as before, we will have to deviate from the 
square-law scaling and increase the switching current by 
nearly two times to 1.05 mA. Thus, we need to pump 
an excess current of 1.05 mA - 0.523 mA = 0.527 mA 
in order to maintain the same switching speed. The cor- 
responding magnetization dynamics without considering 
any thermal fluctuations during the switching is shown 
in Fig. ini where we have clearly recovered the 1.05 ns 
delay. The energy dissipation (dominated by I^Rt) now 
goes up by a factor of two [Is increases by a factor of two 
while T decreases by a factor of two]. Thus, we find that 
if we wish to maintain a constant switching delay, then 
we need to inject some excess current over that dictated 
by square-law scaling and therefore suffer some excess 
energy dissipation. This excess energy dissipation is suf- 
ficiently small so that there is still considerable energy 
saving accruing from the reduction in Ms- Reducing Ms 
by a factor of ~2 results in a net energy saving of ~^14 
times, instead of the ~28 times estimated without impos- 
ing the requirement of constant switching delay. The im- 
portant point is that we have extracted a very significant 
energy saving by reducing Mg by a factor of 2, without 
sacrificing switching speed. 

For illustrative purposes, we show in Fig. [TU] the mag- 
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FIG. 7. Switching dynamics of the magnetization vector in a 
nanomagnet of major axis a — 150 nm, minor axis b = 100 
nm, and Ms = 8 x 10^ A/m, and in-plane shape anisotropy 
energy barrier of 0.8 eV. This simulation does not consider any 
thermal fluctuations during the switching, however, the initial 
orientation of the magnetization is assumed to be Oinit = 
175.5° and (pinit = 90° (thermally mean values for 300 K) 
to avoid the stagnation point exactly along the easy axis. 
Switching is caused by spin-transfer torque induced with an 
in-plane current of 2 mA with 80% spin polarization, (a) polar 
angle 9{t) versus time, and (b) the trajectory traced out by the 
tip of the magnetization vector during switching. Switching 
delay and energy dissipation are 1.05 ns and 1.25 x 10^ kT 
[at room temperature], respectively. 
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FIG. 8. Switching dynamics for Ms = 4.09 x 10^ A/m and 
spin-transfer torque switching current Is of 523 /lA with 80% 
spin polarization. This simulation does not consider any ther- 
mal fluctuations during the switching, however, the initial ori- 
entation of the magnetization is assumed to be Oinu = 175.5° 
and (l>init ~ 90° (thermally mean values for 300 K) to avoid 
the stagnation point exactly along the easy axis, (a) Polar 
angle 9{t) versus time, and (b) the trajectory traced out by 
the tip of the magnetization vector while switching occurs. 
The switching delay and the energy dissipation are 2.1 ns and 
4.3 X W kT [at room temperature], respectively. 



netization dynamics in the presence of thermal fluctua- 
tions at 300 K for the same parameters as in Fig. [HI This 
is one representative run picked out from 10,000 simula- 
tions of the switching trajectory. Note that there is only 
some quantitative difference, but not much qualitative 
difference, between Figs. [^ and [TUl The ripples are some- 
what larger in amplitude and the precessional motion is 
slightly exacerbated. The switching delay has increased 
by ^12% in the presence of thermal agitations, however, 
it should be pointed out that the switching delay may 
decrease as well when the net effect of thermal agitations 
aids the magnetization rotation. 

Fig. [TT] shows how the standard deviations in switch- 
ing delay and energy dissipation due to thermal fluctu- 
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FIG. 9. Switching dynamics for Ms = 4.09 x 10^ A/m and 
spin-transfer torque switching current I^ of 1.05 mA with 80% 
spin polarization. This simulation does not consider any ther- 
mal fluctuations during the switching, however, the initial ori- 
entation of the magnetization is assumed to be Oi„it = 175.5° 
and (pinit = 90° (thermally mean values for 300 K) to avoid 
the stagnation point exactly along the easy axis, (a) Polar 
angle 6(t) versus time, and (b) the trajectory traced out by 
the tip of the magnetization vector while switching occurs. 
The switching delay and the energy dissipation are 1.05 ns 
and 8.6 x 10'' kT [at room temperature], respectively. 
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FIG. 11. Standard deviation in switching delay and energy 
dissipation due to thermal fluctuations at 300 K as a function 
of saturation magnetization Ma in a nanomagnet with fixed 
in-plane shape anisotropy energy barrier of 0.8 eV. (a) Stan- 
dard deviation in switching delay, (b) Standard deviation in 
energy dissipation. 



FIG. 10. Switching dynamics at 300 K for Ms = 4.09 x lO'^ 
A/m and spin-transfer torque switching current Is of 1.05 
mA with 80% spin polarization. The initial orientation of the 
magnetization is 9i„it — 175.5° and (fjimt ~ 90°, which are 
the thermal mean values. This is one specific run from 10,000 
MC simulations, (a) Polar angle 9{t) versus time, and (b) the 
trajectory traced out by the tip of the magnetization vector 
while switching occurs. The switching delay and the energy 
dissipation are 1.18 ns and 9.6 x 10^ kT [at room temperature], 
respectively. 



ations depend on the saturation magnetization Ms- As 
expected, the standard deviation in switching delay in- 
creases with decreasing Ms, because the random ther- 
mal fields hi{t) (i=x,y,z), which are responsible for the 
standard deviation, has a l/^/Mg dependence [see Equa- 
tion ([12])] ■ Furthermore, if we scale Is as M^, then the 
spin-transfer torque also decreases as we reduce Ms and 
that makes the increased thermal field even more effec- 
tive in randomizing the switching delay. For this rea- 
son, the error probability (or switching failure rate) in- 
creases when Ms decreases. This problem too can be 
overcome with some excess switching current. Our sim- 
ulations have shown that if we increase the switching 
current from 523 /xA to 1.05 mA, while holding Ms con- 
stant at 4.09x10^ A/m, then the standard deviation in 
the switching delay goes down from 0.48 ns to 0.23 ns. 
Note that in this way we have recovered approximately 



the same standard deviation in switching delay as that 
for M^ = 8x10^ A/m. 

The standard deviation in energy dissipation however 
shows the opposite trend, i.e. it decreases with decreasing 
Ms- This happens because the energy dissipation Etotai 
is dominated by I^Rt, which is proportional to M*Rt. 
Lowering Ms increases the standard deviation in r, but 
that increase is more than offset by the lower value of 
Ms, so that the net standard deviation in IgRr actually 
decreases with decreasing Ms- The excess current that 
we pump now has a deleterious effect. If we increase the 
switching current from 523 fj,A to 1.05 mA while holding 
Ms constant at 4.09x10'* A/m, then the standard devi- 
ation in energy dissipation goes up from 2.6x10^ kT to 
9.2x10^ kT. 

The ratio of the standard deviation to the mean is 
relatively independent of Ms for both the switching delay 
and the total energy dissipation at 300 K. This ratio does 
not vary by more than 5% when Ms is varied between 
4.09x10^ and 8x10^ A/m. 



IV. DISCUSSIONS AND CONCLUSIONS 

We have shown that one can significantly reduce en- 
ergy dissipation in spin-transfer torque driven switching 
of shape-anisotropic nanomagnets by reducing the sat- 
uration magnetization of the magnet with appropriate 



material choice, while maintaining a constant in-plane 
shape anisotropy energy barrier (by increasing the mag- 
net's aspect ratio) and a constant mean switching speed 
(by pumping some excess current) in the presence of ther- 
mal fluctuations. Also, the increased variance in switch- 
ing delay due to lower saturation magnetization can be 
mitigated by pumping some excess current. In the end, 
by employing these strategies, one can make the energy 
dissipation in spin-transfer torque driven switching of 
nanomagnets competitive with other technologies, with- 
out sacrificing switching speed. This bodes well for ap- 
plications of spin-transfer torque switched nanomagnets 
in non-volatile logic and memory. 
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O ; SI Fluctuation about the easy axis due to thermal torque 

l> 

^^ i We here explicitly show that thermal fluctuations by themselves can dislodge the magnetization 
from the easy axis. At T = K, 9init = 180°. In this case, spin-transfer-torque can never budge 
the magnetization vector since it vanishes whenever sin 6 = 0. However, when T > K, thermal 
fluctuations can dislodge the magnetization from the easy axis. 

K> ' In the spherical coordinate system, with constant magnitude of magnetization 



S 



Accordingly, 



^^ = ^^e, + .zn^(t)^^e,. (SI) 

a ( n„i(i) X "^"""^^M = -asine{t) (l)'{t) e^ + aO'it) e^ (S2) 



and 

dn^{t) ( dn^{t) 



dt 



a (n^{t) X !^^}b^\ = [e'{t) + asin9{t) (t>'{t)] ee + [sin9{t) (j)'{t) - aO'it)] e^ (S3) 
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where ()' denotes d{)/dt. Equating the eg and e^ components of the LLG equation [Equation (16) 
in the main paper] and using Equations (8) and (10) in the main paper, we get 

9'{t) + asine{t) (p'it) = -^[Soe(t) sine{t) -ssine{t) + PJt)] (S4) 

My 

sine{t)(j)'{t)-ae'{t) = ^[2B{t)sine{t)cose{t) + s sine{t) - Pe{t)]. (S5) 

My 

When sin 6 = 0, Equations (|S4p and (jSSp yield 

0'{t) = ^P^t) (S6) 

aO'it) = ^Pe{t). (S7) 

Substituting for Pei^t) and P(j,{t) from Equations (14) and (15) as in the main paper and using 
d = 180°, we get 

ahx{t)sin(l){t) — ahy{t)cos(t){t) = hx{t)cos(p{t) + hy{t)sin(t){t) (S8) 

which yields 

Using this value of (\>{t) in either Equation (|S6|) or Equation ()S7p . we get 

'^{ahy{t) + K{t)Y + (a/i,(t) - hy{t)Y 

We can see from the above equation that thermal torque can dislodge the magnetization vector 
from the easy axis because the time rate of change of 9{t), i.e. 0'{t), is non-zero, even though 
sinO = 0. Note that the initial rate of deflection of the magnetization vector from the easy axis 
does not depend on the component of the random magnetic field along the z-axis [/i2(t)], which is 
a consequence of having the z-axis as the easy axis of the nanomagnet. However, as soon as the 
magnetization is deflected from the easy axis {sinO ^ 0), all the three components of the random 
field would come into play. 

S2 Thermal distributions for Oinu and (jyinit 

Fig. IS II shows the fluctuations in the polar angle 9{t) over time due to thermal agitations at 300 
K, when the magnetization is nominally along one of two stable orientations along the easy axis {0 
= 180°). There is no external torque (e.g., spin-transfer torque) here, but there is a torque due to 
shape anisotropy since it is a torque internal to the nanomagnet. 

In Fig. ISl) there are 1 million time steps in the 100 ns interval and this time interval is long 
enough to make the mean value of 6{t) independent of the interval duration. Fig. [S2] shows the 
distributions of the polar angle 0{t) and the azimuthal angle (j){t) due to thermal fluctuations at 
300 K within the same time interval of 100 ns used in Fig. ISII The most likely value of is 180° 
which is a stagnation point where spin-transfer-torque vanishes, but the mean value is not 180°. 
The mean value is ~175.5°, so that the time-averaged (mean) deviation from the easy axis is -^4.5°. 
We noticed that halving the value of Mg does not have any effect on the mean value oi 9, or the 
deviation from the easy axis, because the in-plane shape anisotropy energy barrier is kept constant 
by adjusting the three components of demagnetization factor. Even though variance of the thermal 
field h{t) depends on M<j, its mean value is zero, and hence the mean value of the deviation from 
the easy axis turns out to be relatively independent of Mg. The distribution of 9{t) appears to be 
nearly exponential, while the distribution of (j){t) contains two nearly Gaussian distributions peaked 
at the two in-plane angles (j) = ±90°. Unlike the ^-distribution, the (/(-distributions are symmetric 
about the mean values. 
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Figure SI: Temporal fluctuations in the polar angle 6{t) around one of the easy axes [6 = 180°) 
due to thermal torque at room temperature. The torque due to shape anisotropy is present since 
it is internal, but there is no spin-transfer-torque which is always applied from an external source. 
Note that since 6 is the polar angle, it is constrained to the interval [0°, 180°]. The nanomagnet 
has saturation magnetization Mg = 8 x 10^ A/m, the major axis o = 150 nm and the minor axis b 
= 100 nm. The in-plane shape anisotropy energy barrier is 0.8 eV or ~32 kT at room temperature. 
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Figure S2: Distributions of the polar angle 0{t) and azimuthal angle (j){t) due to thermal fluctuations 
at room temperature in a nanomagnet of elliptical cross-section. The major axis a = 150 nm, the 
minor axis h = 100 nm and the saturation magnetization Mg = 8 x 10^ A/m. The in-plane shape 
anisotropy energy barrier is 0.8 eV or ~32 kT at room temperature. The magnetization vector is 
assumed to point in one of two stable orientations along the easy axis {9 = 180°) but is perturbed 
by thermal fluctuations, (a) Distribution of polar angle 9[t) due to thermal fluctuations at room 
temperature. The mean of the distribution is '^175.5°. (b) Distribution of azimuthal angle (t){t). 
There are two peaks in the distribution centered at 90° and 270° (or -90°) that correspond to the 
plane of the nanomagnet. 
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Distributions of switching delay and energy dissipation in the 
presence of thermal fluctuations 



Fig. IS3l shows the room temperature (300 K) distributions of switching delay and energy dissipation, 
respectively, when both spin-transfer torque and thermal torque act on a nanomagnet with in-plane 
shape anisotropy energy barrier of 0.8 eV. The spin-transfer torque is due to an in-plane current 
of 523 /xA with 80% spin polarization. The saturation magnetization Ms = 4.09 x 10^ A/m. Here, 
10,000 trajectories were simulated to compute these distributions. We assume that when a spin- 
polarized current is applied to initiate switching, the magnetization vector starts out from near the 
south pole {0 ~ 180°) with a certain {9init,4>init) picked from the initial angle distributions at 300 
K. Switching is deemed to have completed when 0[t) becomes equal to or less than 4.5°. Note that 
both distributions are asymmetric about the mean. Regrettably, the long-delay tail extends much 
farther than the short-delay tail. The high-delay tail is associated with those switching trajectories 
that start very close to ^ = 180° which is a stagnation point. In such trajectories, the starting 
torque is vanishingly small, which makes the switching sluggish at the beginning. During this time, 
switching also becomes susceptible to backtracking because of thermal fluctuations, which increases 
the delay further. The distribution of energy dissipation is complex, which shows multiple peaks 
in the distribution. The distribution is due to the variation in internal energy dissipation, E^ (see 
Fig. 6 in the main paper). 
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Figure S3: Room temperature (300 K) distributions of switching delay and energy dissipation when 
a spin-transfer-torque is applied to switch a nanomagnet from one mean thermal deviation from one 
orientation along the easy axis to one mean thermal deviation from the other diametrically opposite. 
Here, the saturation magnetization Mg = 4.09 x 10^ A/m and the in-plane shape anisotropy energy 
barrier is 0.8 eV. The major axis of the nanomagnet is a = 300 nm and the minor axis is 6 = 50 
nm to cause this shape anisotropy energy barrier. The magnitude of the in-plane spin-polarized 
current generating spin-transfer-torque is 523 /xA and the spin polarization of the current is 80%. 
(a) Distribution of switching delay: mean value = 2.1 ns, standard deviation = 0.48 ns, and (b) 
distribution of energy dissipation: mean value = 4.1x10^ kT [at room temperature], standard 
deviation = 2.6x10^ kT [at room temperature]. 



S4 The beneflts of switching with a larger switching current /<, 

It is obvious that increasing switching current (and hence energy dissipation), while keeping satu- 
ration magnetization and shape anisotropy energy barrier constant, will result in faster switching. 
In order to illustrate this point, we do not consider any thermal fluctuations during the switch- 
ing to avoid complications; however, the initial orientation of the magnetization is assumed to be 
Oinit = 175.5° and (t>init = 90° to avoid the stagnation point exactly along the easy axis. Fig. [S4] 
shows the switching dynamics for this case when switching is caused by spin-transfer-torque induced 
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with an in-plane switching current Ig of 2 mA possessing 80% spin polarization. The saturation 
magnetization Ms = 8 x 10^ A/m. The switching occurs in ~1 ns, while dissipating 1.25x10^ kT 
[r = 300 K] of energy. If we increase Is to 4 mA (Fig. ISSP while keeping everything else the same, 
the switching delay r drops to 455 ps, but of course the energy dissipation (dominated by I'^Rt) 
increases by 75%. What Figs. IS4l and IS5I show is that the reason a higher switching current decreases 
switching delay is because it suppresses the ripples in the transient dynamics of the magnetization 
vector. The ripples are caused by complicated precessional motion of the magnetization vector 
shown in Fig. 94(b)] A larger Is exerts a larger spin-transfer torque on the magnetization vector 
that suppresses its errant precessional motion and brings about the switching faster. The impor- 
tant point is that since one can always decrease switching current (and hence energy dissipation) 
by sacrificing switching speed, it is imperative to keep the thermal mean of the switching delay (r) 
constant when studying how the switching current Is scales with saturation magnetization Mg. 



200 r 



150 



Dynamics of e 



Dynamics of magnetization axis 



S?100 




0.6 Of 
Time (ns) 

(a) 



1.2 



14 




-1 -0.1 



(b) 



Figure S4: Switching dynamics of the magnetization vector in a nanomagnet of major axis a = 150 
nm, minor axis h = 100 nm, and Mg = 8 x 10^ A/m, and in-plane shape anisotropy energy barrier of 
0.8 eV. This simulation does not consider any thermal fluctuations during the switching, however, 
the initial orientation of the magnetization is assumed to be Oinit = 175.5° and (jjinit = 90° (thermally 
mean values for 300 K) to avoid the stagnation point exactly along the easy axis. Switching is caused 
by spin-transfer torque induced with an in-plane current of 2 mA with 80% spin polarization, (a) 
polar angle 9{t) versus time, and (b) the trajectory traced out by the tip of the magnetization vector 
during switching. Switching delay and energy dissipation are 1.05 ns and 1.25 x 10^ kT [at room 
temperature], respectively. 



Increasing the switching current, while keeping everything else constant, also decreases the 
spread (variance) in switching delay caused by thermal fluctuations at non-zero temperatures. 
Monte Carlo simulations show that the standard deviation in switching delay decreases from 0.23 
ns to 0.14 ps at room-temperature when the switching current is increased from 2 mA to 4 mA. 
This happens because the larger switching current suppresses precessional motion of the magneti- 
zation vector that subjects the switching dynamics to greater variability in the presence of thermal 
fluctuations. A smaller standard deviation in the switching delay helps increasing the operational 
clock-frequency and facilitates clock synchronization on a chip. Thus, expending more energy to 
switch (larger Is) bears some distinct advantages in memory and logic applications, but of course 
at the cost of increased energy dissipation. 

In Fig. IS61 we plot the dynamic switching error probability as a function of switching delay 
r for different switching currents Ig in the presence of thermal fluctuations at room temperature. 
Switching error occurs when the magnetization vector sets out from its initial stable orientation 
along the easy axis towards the intended stable orientation diametrically opposite, but fails to reach 
the latter within the stipulated delay r because of thermal agitations. For a fixed switching delay, 
a higher current ensures a higher probability of successful switching because it provides a stronger 
spin-transfer-torque that overwhelms the thermal torque and suppresses the errant precessional 
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Figure S5: Switching dynamics of the magnetization vector in a nanomagnet of major axis a = 
150 nm, minor axis b = 100 nm, and Mg = 8 x 10^ A/m. The resulting in-plane shape anisotropy 
energy barrier is 0.8 eV. This simulation does not consider any thermal fluctuations during the 
switching, however, the initial orientation of the magnetization is assumed to be Oinit = 175.5° and 
(t'init = 90° (thermally mean values for 300 K) to avoid the stagnation point exactly along the easy 
axis. Switching is caused by spin-transfer torque induced with an in-plane current of 4 mA with 
80% spin polarization, (a) polar angle 9{t) versus time, and (b) the trajectory traced out by the 
tip of the magnetization vector during switching. Switching delay and energy dissipation are 455 
ps and 2.2 x lO'^ kT [at room temperature], respectively. 
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Figure S6: Dynamic switching error probability as a function of switching delay for different switch- 
ing currents at T = 300 K. The nanomagnet has dimensions of o = 150 nm, h = 100 nm, and the 
saturation magnetization is Mg = 8 x 10^ A/m. The resulting in-plane shape anisotropy energy 
barrier is once again 0.8 eV. Faster switching with a constant current increases the likelihood of 
error (failure to switch), but at a fixed speed of switching, the error probability decreases with 
increasing switching current. 
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motion of the magnetization vector that increases the chances of backtracking and the resulting 
switching failure. For a fixed switching current, the error probability increases with the speed of 
switching. This is somewhat obvious; if we allot less time for the switching to complete, it is less 
likely that switching will complete in that allotted time. 

S5 Switching delay-energy as a function of Ms 

Fig. IS7I shows the switching delay-energy trade-off as a function of saturation magnetization Mg . 
The plot shows a clear trade-off between the switching delay and energy dissipation as Ms is varied. 
As Ms is lowered, less energy is dissipated, however, switching delay is increased. The square law 
scaling of switching current Ig oc M^ is used. As discussed in the main paper, with lower values 
of Ms, a higher magnitude of current (than that of using the square law Is oc Mg) can reduce the 
switching delay to some smaller value corresponding to a higher value of Ms , while dissipating lesser 
energy than that of the high- valued Ms . 
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Figure S7: Thermal mean of the total energy dissipation versus thermal mean of the switching delay 
as a function of saturation magnetization Ms (4.09 x 10^ — 8 x 10^ A/m). The square law scaling of 
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is increased. 



switching current Is cc Ms is used. A lower Mg incurs less energy dissipation but switching delay 



